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Abstract Two-particle femtoscopy reveals the space-time substructure of the 
freeze-out configuration from heavy ion collisions. Detailed fingerprints of bulk col- 
lectivity are evident in space-momentum correlations, which have been systemati- 
cally measured as a function of particle type, three-momentum, and collision condi- 
tions. A clear scenario, dominated by hydrodynamic-type flow emerges. Reproduc- 
ing the strength and features of the femtoscopic signals in models involves important 
physical quantities like the Equation of State, as well as less fundamental techni- 
cal details. An interesting approximate "factorization" in the measured systematics 
suggests that the overall physical freeze-out scale is set by final state chemistry, but 
the kinematic substructure is largely universal. Referring to previous results from 
hadron and lepton collisions, we point to the importance of determining whether 
these "universal" trends persist from the largest to the smallest systems. We review 
theoretical expectations for heavy ion femtoscopy at the LHC, and point to direc- 
tions needing further theory and experimental work at RHIC and the LHC. 
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Amongst all the classes of observables analyzed in subatomic collisions, fem- 
toscopy is unique as it directly addresses the space-time structure of the evolving 
system at the Fermi scale. It exploits the fact that the likelihood to emit particle 
a with momentum is modified by the emission of particle b with momentum 
Ph, so long as there is any structure- due to quantum (anti-)symmetrization and/or 
final state interactions- in the two-particle relative wavefunction ^{q^r), where 
q = [pa — pi,)/2 and r is the relative space-time separation. The (^-dependence of 
the modification, which can be measured, reveals the separation distribution of the 
particles. 

The techqnique of measuring the arrival-time correlation of radio-wavelength 
photons to extract angular sizes of stars was invented more than 50 years ago by 
Hanbury Brown and Twiss (US)- Although the entire approach was initially greeted 
skepticall}0, the technique was used for years in the ground-based astronomy com- 
munity. 

Five years later, correlations of identical pions were shown to be sensitive to 
source dimensions in proton-antiproton collisions by Goldhaber, Goldhaber, Lee, 
and Pais (GGLP) (jsl). In the 1970s, these methods were refined by Kopylov and 
Podgoretsky (0; □), Koonin, (@), and Gyulassy 

i), 

and other classes of correlations 
were shown to be useful for source-size measurements, such as strong and Coulomb 
interactions. 

While the "GGLP effect" in particle physics and the original "HBT effect" are 
in principle based on different theoretical frameworks (13 U), in practice the ex- 
tracted information and techniques are quite similar. Thus, when in the 1980's ex- 
periments at the Bevalac (Il2tll3h made clear that identical pion correlations yielded 
quantitative space-time information in nuclear collisions, the heavy ion community 
referred to such measurements as "HBT analyses." Around the same time, the tech- 
nique began to be applied to prot on-proton and fragment-fragment correlations in 
low energy nuclear collisions \ 

HBT measurements were among the first results reported when studies of truly 
relativistic heavy ion collisions began twenty years ago with the availability of Oxy- 
gen and Silicon beams SPS/CERN and AGS/BNL, respectively. Since then, extract- 
ing space-time scales from two-particle correlations has become an integral part 
of the heavy ion program worldwide. The availability of high-statistics data from 
large-acceptance detectors and an continuing development of theoretical formalism 
has tremendously expanded the range of geometrical information extracted from the 
correlations. Moving far beyond the original "source size" measurements, modern 
analyses routinely extract the extent, shape, orientation, and dynamical timescale 
of the emission zone. This information has been systematically extracted and com- 



^ Towards the end of his Hfe, Hanbury Brown gave a retrospective presentation {J on this seminal 
work. In its entirety, his listed abstract read "The talk will give a brief history of the early develop- 
ment of Intensity Interferometry and its subsequent battle against common sense." See also a. 
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pared as a function of collision energy, centrality and ion species. It is further binned 
in terms of pj, y and (/), and for a growing variety of particle species. 

Data systematics and model comparisons to RHIC, SPS, and AGS data have re- 
cently been reviewed in detail Referring the reader there for full details, we 
wish here to present a minimally self-contained overview of methodology and sta- 
tus of femtoscopy in heavy ion collisions on the eve of the LHC. We will discuss 
generic features of the source in these collisions, which we consider firmly estab- 
lished. Finally, we will briefly discuss some expectations for femtoscopy at the LHC, 
and open questions to be addressed there. 



1.1 The urgency of space-time in heavy ion studies 

The goal of the relativistic heavy ion physics program is the creation and study 
of color-deconfined bulk matter. The hope is to to better understand and test QCD 
in the sector in which it is 1) most interesting and at the same time 2) most diffi- 
cult to calculate. It is most interesting in the low momentum transfer (alternatively 
large length scale or "soft") sector, since it is here where its most characteristic and 
important feature- color confinement- is manifest. On the other hand, it is most 
theoretically difficult since the coupling between the colored objects in the theory 
is strong, rendering perturbative techniques ineffective. This feature has a major up- 
side, however. In particular, it gives hope that heavy ion collisions might produce a 
bulk system, at least partially thermalized, which can be described in the language 
of bulk quantities like temperature and pressure. 

The most celebrated early result from the RHIC program was the success of hy- 
drodynamic calculations to describe the transverse momentum and azimuthal struc- 
ture of measured hadron spectra. While the focus has since shifted to more refined 
and important details related to initial conditions, pre-equilibrium collectivity, fi- 
nite viscosity and equations of state, there remains near-universal consensus that 
the measured soft-sector momentum distributions are clear evidence of bulk collec- 
tive behaviour- hydrodynamic flow. However, projection of the seven-dimensional 
(space-time and momentum) emission pattern onto momentum-only axes represents 
only a minimal test of the complex, non-trivial geometric substructure predicted by 
hydrodynamics. In particular, collective flow is much more precisely characterized 
by signature correlations between momentum and coordinate space. Femtoscopy 
provides the most direct experimental access to such patterns. A solid case for a 
bulk system in heavy ion collisions depends crucially on such measurements. 

Intrinsic thermodynamic quantities which characterize a state of matter, require 
knowing the volume of the sample of matter under study. The relationship between 
these quantities constitutes the Equation of State (EoS). Although much of the infor- 
mation we seek is for the static, equilibrium case, experimental information on QCD 
in this sector can only be obtained through the study of highly dynamic systems- 
i.e. collisions. Furthermore, the dynamics reflect the equation of state and viscosity. 
Without the spatial information, inferences of bulk quantities like pressure, energy 
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density or entropy density would require untested assumptions about the volume, 
and without the temporal information, inferences into dynamics would be largely 
faith-based. 

The high inherent dimensionality to correlation functions (six dimensions for 
each species pair and for each beam energy, beam/target combination and central- 
ity) make correlations the ultimate test for any model that claims to describe the 
evolution of a heavy ion collision. Thus, it is no surprise that femtoscopic data have 
proven the most challenging to explain of all heavy-ion observables. For these rea- 
sons femtoscopy is perhaps the most important soft-sector observable in relativistic 
heavy ion physics. 



1.2 The neccessity of systematics in heavy ion studies 



The physics program at the Relativistic Heavy Ion Collider (RHIC) at Brookhaven 
National Laboratory is remarkably rich, thanks to the machine's unique ability to 
collide nuclei from to '^^Au, in fully symmetric (e.g. Ah-A or pH-p) to strongly 
asymmetric (e.g. d+Au) entrance channels, over an energy range spanning more 
than an order of magnitude. The capability to collide polarized protons provides 
access to an entirely new set of fundamental physics, not discussed further here. 

Achieving the primary aim of RHIC- the creation and characterization of a color- 
deconfined state of matter and its transition back to the confined (hadronic) state- 
requires the full capabilities of RHIC. In particular, comparisons of particle distri- 
butions at high transverse momentum {pj) from Ah-A and pH-p collisions, probe the 
color-opaque nature of the hot system formed in the collisions ([IlliTlB. Com- 
parison with reference d + A collisions were necessary to identify the role of initial- 
state effects in the spectra (Il9h . Comparing anisotropic collective motion from non- 
central collisions of different-mass initial states (e.g. AuH-Au versus Cu + Cu) (l20l) 
tests the validity of transport calculations crucial to claims of the creation of a "per- 
fect liquid" at RHIC (|2I[) . Indeed, a main component of the future heavy ion pro- 
gram at RHIC involves a detailed energy scan, designed to identify a predicted crit- 
ical point in the Equation of State of QCD (1221) . 

The need for such systematic comparisons is not unique to RHIC, but has been 
a generic feature of all heavy ion programs (l23t l24r) . from low-energy facilities like 
the NSCL (Michigan State), to progressively higher-energy facilities at SIS (GSI), 
the Bevati-on/Bevalac (Berkeley Lab), AGS (Brookhaven), and SPS (CERN). The 
nature of heavy ion physics is such that little is learned through study of a single 
system. 

The femtoscopic excitation function, in particular, was expected to exhibit an un- 
ambiguous signal of a hoped-for sharp phase transition, based on general grounds 
and essentially independent of a model ( 1251 : l26c 1271) . The absence of such a signal 
was originally termed the "RHIC HBT puzzle" (l28h . However, our understanding 
of the deconfinement transition at RHIC has evolved away from sharp first-order 
phase transition, thanks in part to the femtoscopic measurement. In addition, more 
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detailed refinements in the comparison of models with measured correlation func- 
tions, discussed below, have largely ameliorated this puzzle. Nevertheless, we will 
mention further possible "puzzles" in direct comparisons of femtoscopy from heavy 
ion versus elementary particle collisions. 



1.3 Organization of this paper 

In Section |2] we present a brief primer on two-particle femtoscopy, sketching the 
theoretical formalism and connection of the correlation function to the structure of 
the source function. In Section[3] we review the systematics mapped out by a broad 
range of femtoscopic measurements, and the scenario which clearly emerges. There, 
we also briefly discuss what has been observed in high energy hadron-hadron and 
lepton-induced collisions. In Section |4] we discuss conclusions drawn from para- 
metric and transport model comparisons so far. In Section|5] we review theoretical 
expectations for heavy ion collisions at the LHC, expected in just a few years. In 
Section|6] we summarize. 



2 A Femtoscopy Primer 

Here, we give a very brief overview of the formalism of femtoscopic correlation 
functions. For a more complete discussion of the theoretical and experimental de- 



tails, see ( 15). 



2.1 What is measured - correlation functions and source functions 

Femtoscopic analyses are based on two-particle correlation measurements at small 
relative momentum. For particles of types a and b, the six-dimensional correlation 
function is typically expressed in terms of the relative and total momentum, 

^^^''^^ {dN"/d^p^,){dN>'/d^p,)' 

P^Pa + Pt, q^^{pa-pbr/2- ^''"'J^^'^ Pr 

The last term in the definition of the relative momentum q projects out the temporal 
component in the two-particle rest frame. At times, we will use Q = 2qto refer to 
the relative momentum. 

Experimentally, the formal definition of the correlation function in Equation[T]is 
seldom used in heavy ion physics. Instead, it is approximated by 
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where Ap'(q) is the signal distribution, Bp'{q) is the reference or background distri- 
bution which is ideally similar to A in all respects except for the presence of femto- 
scopic correlations. In heavy ion collisions, B is usually formed via the event-mixing 
technique (29), in which q (or any pair-wise variable) is formed when particles a and 
b come from different events; in contrast, the histogram Ap'(q) represents the dis- 
tribution when the particles come from the same event. Except in measurements in 
which the acceptance is extremely limited (13) or for very low multiplicities Jsol) . 
Equations [U and [2] yield practically identical results IbTI). In any event, the refer- 
ence distribution B should account for effects of detector acceptance and efficiency, 
and fluctuaions in reaction-plane orientation, multiplicity, and vertex position; ac- 
counting for such effects are most naturally done with the event-mixing technique. 
(^p(q) represents a potential correction factor introduced to compensate for non- 
femtoscopic correlations present in the signal that are not fully accounted for in the 
background as well as artifacts resulting, e.g., from finite resolution and contamina- 
tion. See (iTsI) and references therein, for a more detailed discussion. 

The designation of the total momentum P by a subscript in Equation |2] reflects 
the fact that usually experimental correlation functions are formed in terms of q and 
analyzed separately for various selections of the total momentum P. As discussed in 
Section[Tl femtoscopic correlations arise only through the spatially-dependent two- 
particle relative wavefunction 0„ f,(q,r), which vanishes as |q oo\. Thus, for large 
relative momentum, the correlation function approaches a constant value; typically, 
C(|q ^ oo|) is normalized arbitrarily to unity. The presence of other correlations 
(e.g. due to momentum conservation dsj) ) can complicate matters, especially for 
small systems. 

Femtoscopic correlations are expected to be described by the Koonin equation 
i), 

Jd'r' S"'{Py)\0{q'y)\\ (3) 

Here, the squared relative outgoing two-particle wave function serves as a weight 
which enhances or reduces the probability for emitting two particles depending on 
their separation. In Eq. (O r' refers to the separation of the two particles in the two- 
particle rest frame and q' is the relative momentum in that frame. The probability of 
emitting a single particle of momentum p from space-time point x is s{p,x), and is 
evaluated at a momentum such that each particle has the same velocity and such that 
the total momentum is P. Justification of the Koonin equation requires assuming that 
the sources are highly uncorrected (chaotic), that the phase space density is not too 



C''^P,q) = 

cab Co U\ = 
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large, and that the relative momentum is small. A detailed discussion is provided in 
(El. 

Since the separation of two particles of the same velocity is independent of time, 
the source function can also be written in terms of the asymptotic phase space den- 
sities in the pair frame (where both momenta are zero), 

s.h(p ^ Jd'^'gdhj u (p;. - o,x;„OA (p; = o,x;,o 5' (r- - [x^, - x;]) 
/«'V„flf3^;,/„(p^, = o,x;„f)A(p; = o,x;,f) 

assuming t is any time after all particles have been emitted. It is more accurate to 
state that femtoscopy measures the structure of outgoing phase space clouds of a 
given velocity, rather than claiming that the source is measured. 

Since correlations are six-dimensional objects, there is no hope to extract all 
seven dimensions of information in i(p,x). The separation x' = 7(x — vpf ), the spa- 
tial separation along the direction of P, is inextricably linked to the time. Nonethe- 
less, if one assumes that the spatial separation along P is distributed similarly to 
one of the transverse separations, one can infer the temporal information by con- 
sidering the difference. Such assumptions usually become trustworthy only for long 
lifetimes, as has been measured femtoscopically in the case for heavy ion collisions 
at low energy S m S), < 50A MeV. We discuss this further in Section |3] in 
relation to the ratio Rout /Rside- 

Furthermore, since one can only infer information about the relative separation 
of two particles with a given velocity, one measures a subset of the entire source, 
even if the correlation measurements span all of momentum space. The subset of 
the volume relevant for femtoscopy is referred to as the region of homogeneity ( l36t 



151) . and can strongly vary with P, especially in the presence of collective flow, as 
discussed in Section[3] 

From viewing Equations[3]and|4]it is clear that correlations can, at best, provide 
the source function S{P,r), and the principal role of femtoscopy is to discern all 
possible information about the source function from the correlation function. 

Rather than expressing the correlation function in terms of the three components 
of q or the source function in terms of the three components of r, one can decompose 
the ang ular information using a basis of either Cartesian or spherical harmonics (l37t 



38tl32h . In Cartesian harmonics the components are labeled by ^ = {£x,iy,iy), with 
i = £x + £y + iz^ ^nd the decomposition of an angular function F{Q) is defined by 

where Q is a unit vector Some examples of the Cartesian harmonics, A({Q), are: 

Ao.0,0 = lAiM^) = axAififii^) - - (l/3),Ai,Ko(i^) = Av4.,---. The 
convenience of such decompositions derives from the property that a given angu- 
lar component of the correlation function has a one-to-one link to the corresponding 
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element of the source function, 

Ct{q) = 1 4nr^drK({q,r)Se{r), (6) 



Kei^^r) = ^ / ^^cos0^r^f(cos0,,,.)|0(^,r,0^,.)l^ 



1 
2 

where O^r is the angle between q and r. Equation (|6]l makes it possible to evaluate 
the correlation coefficients Q component-by-component. Decompositions are also 
an efficient way to view features of the correlation arising from a specific symmetry. 
For instance, in a correlation function of non-identical particles, it is insightful to 
analyze the (^out = l,^side = ^long = 0) projections. Spherical harmonics have the 



same convenient properties, and several simplifying symmetries (1321) 

Any transport theory can generate emission probabilities, and can therefore pro- 
vide a source function from which to calculate correlation functions; details of the 
calculation are discussed in (jTsl). Comparisons of the calculated correlation func- 
tions to the measured ones represent the best possible test of the space-time struc- 
ture of the model's freezeout configuration. However, it is often useful (for example 
when the calculated and measured correlation functions disagree) to parameterize 
source function to identify its dominant features, such as resonance tails, long life- 
times, or anisotropic shapes. The extracted parameters, which we discuss next, are 
typically what are published and compared between experiments or to theory. 



2.2 Representation of the source function 

Here, we briefly review parameters and formalism commonly used to present and 
summarize geometric information from femtoscopic measurements. 



2.2.1 Reference frames and coordinate systems 

The most common coordinate system in which to decompose the relative momen- 
tum is the Bertsch-Pratt system (i39i : i4Q) . in which q ~ {qoutiqsideiliong}- The "long" 
axis is oriented along the beam direction; the transverse momentum vector Pj (i.e. 
the projection of the total momentum perpendicular to the beam) defines the "out" 
direction; the "side" direction is perpendicular to these. These are indicated in Fig- 
ure H 

The correlation function is generally presented in one of three reference frames. 
The Longitudinally Co-Moving System (LCMS) is defined such that the pair's total 
momentum along the beam direction vanishes (P^ ~ 0). This is commonly used at 
ultra-relativistic energies, for which the system is approximately boost-invariant. In 
a boost-invariant scenario, pairs with a given velocity in the beam direction V; are, on 
average, emitted from source elements moving with the same V,. Since comparisons 
to e+ — collisions are increasingly interesting, as we discuss later, one should be 



10 



M.A. Lisa and S. Pratt 



aware that in these measurements, the LCMS system is defined such that the pair 
momentum along the event's thrust axis vanishes (14 lb . 

In some analyses, an additional transverse boost is made to the Pair Center of 
Mass System (PCMS), in which the total momentum of the pair vanishes. This dif- 
fers from the LCMS only by a Lorentz factor 7 in the outward direction, but some 
prefer the PCMS since the relative wave function is most conveniently expressed in 
this frame. For instance, a sharp resonance peak is no longer sharp if the correla- 
tion is viewed away from the pair frame. The PCMS is usually used in analyses of 
non-identical particle correlations (in which the relative momentum is sometimes 
denoted k*) and when performing "imaging" fits as described below. 

At low collision energies, for which boost invariance is a poor approximation, the 
CMS of the colliding system (the lab frame, for symmetric collisions at a collider) is 
often used (e.s.Ali): this would be the natural frame for a source with no longitudinal 
flow. 



2.2.2 Gaussian parameterizations 

To extract the dominant scales and shape of the homogeneity region, a Gaussian 
source function is often used. In general, the source function of the separation r' = 



S(Rr') 



• exp 



- L L ('•;-^0(o-^.)-[(^'rl-,,/4 

i=o,s,l j=o,sJ 



(7) 



where o,5,Z stand for the "out," "side" and "long" vector components. Nine quan- 
tities characterize the distribution: the vector difference between average emission 
points A and a symmetric 3x3 covariance matrix 



(8) 
(9) 



where (• • •) indicates an average over the source function S"''. An implicit depen- 
dence of A and on P has been suppressed. In most analyses, several of these nine 
factors vanish due to symmetry (e.g. l43tl44l) . 
Identical particles 

For identical particles {a = b), A —0 trivially, and the covariance matrix is typi- 
cally expressed in terms of the six "HBT radii parameters" Rjj. For the commonly- 
measured two-pion correlations, they are usually extracted by fitting the measured 
correlation function with the functional form 



Cp (q) = [a/TcouI (Q) f 1 + + (1 - A) 



(10) 



where is the overall normalization discussed above. Here, AcouK^) is a factor, 
which approximately accounts for the Coulomb interaction between pions (45; 46l) . 



Feratoscopically Probing the Freeze-out Configuration 



11 



For fully chaotic emission from a purely Gaussian source, and with perfect exper- 
imental particle identification (PID), the factor A should be unity. If experimental 
identification is perfect and the source is purely chaotic, X represents the fraction 
of correlated pairs, i.e., the fraction of pairs where neither pion is emitted from far 
outside the source being described by the Gaussian. For instance, if 25% of the pi- 
ons were either emitted from long-lived resonances or were mis-identified electrons, 
one would expect the lambda parameter to be 0.75^. Typically A ~ 0.3 — 0.7, and in 
a typical fit to a measured correlation function, the covariance between the different 
"HBT radii" is small or zero, while the covariance between radii and X is usually 
slightly positive. 

The diagonal components R^, and Rj characterize the lengths of the major 
axes of the homogeneity region of particle a, as shown in Figure|3] The off-diagonal 
components R^^, R^/ and R^i quantify the orientation of these axes relative to the 
"out", "side" and "long" directions; they encode important information about fem- 
toscopic anistropies of the collision. The azimuthal dependence of R~^, relative to 



the reaction-plane, probes the transversely anisotropic shape of the source (1471 : 1481) : 
it is the coordinate-space analog of the momentum anisotropy V2- Finally, the az- 
imuthal dependence of probe the spatial "tilt" of the collision zone ( l48h : it is 
the coordinate-space analog of directed flow V]. 



As a dynamic probe, the rapidity dependence of R^^i reveals the extent to which 



the dynamics of the colliding system is boost-invariant (l49h . A finite value of this 
parameter, as measured in LCMS indicates a breakdown of boost invariance. For 
central collisions, an alternate Gaussian source parameterization is sometimes used, 
instead, to probe boost-invariance. Using the Yano-Koonin parameterization 1^ 



5 lb . the correlation function is fit with the form (52) 



Cp{Q)=n\XKcou\ (q) (l +e-''i^i-''H(2r25)-(«o-'^H)(2-")') + (1 - A) 



(11) 

Here, u = 7(1,0,0, vy^:), and vyk is the velocity of the source element, along the 
beam direction, which emitted the pions. In a boost-invariant system this velocity is 
identical to the longitudinal velocity of the pion pair 
Non-identical particles 

Distinguishable particles (a ^ b) with the same velocity can have different av- 
erage space-time freezeout positions, as quantified by the displacement vector A in 
Equation |7] For boost-invariant sources, 4/ vanishes, as does As for analyses inte- 
grated over azimuth. Few measurements of the displacements have been performed 
to date: only displacements along the pair velocity have so far been published at 



low (|53|:|54|:|55|:|56|) and RHIC W3 energies 



2.2.3 Non-Gaussian parameterizations 

While Gaussian scale and covariance parameters capture the bulk properties of the 
emitting source, there is no reason to expect S"'' (r') to be purely Gaussian. Indeed, 
discrepancies have long been noted between measured correlation functions and 
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fits using Equation [T] These discrepancies are not large, but are to be expected. 
For example, exponentially decaying resonance feed-down contributions will natu- 
rally lead to long-range tails in the emission probability, creating a "core-halo" iSab 
structure. The "halo" affects the low-|q| structure of the correlation function. Un- 
til recently, such structure could not be resolved, and the halo mostly affected the 
A factor in Gaussian fits (js^)- However, high-resolution data at RHIC has allowed 
investigation into the source halo. 

Various non-Gaussian functional forms and expansions have been proposed ((SS 
which fit the data better than a pure Gaussian, for example Edgeworth expansions 
(fed) , or double-Gaussian fits. However recently there has been considerable activity 
in "imaging" the source ( 1611 : l62t l63l : l64t l65l: 1661) . As with the Gaussian fits discussed 
above, imaging is essentially a multi-parameter fit of a calculated correlation func- 
tion to the measured one. The calculated correlation function is generated through 
Equation [3] where the source function S"'^{r') is parameterized as a sum of spline 
functions. The weights in this sum are the fit parameters varied to minimize x^- 
Considerable mathematical and numerical machinery has been developed in order 
to achieve a stable fit- including setting bin sizes and stationary points (or "knots") 
by hand- and in practice the technique is nontrivial. It may be fairly said that prop- 
erly imaging a source remains something of an art at this point. Mutli-dimensional 
imaging makes use of the harmonic decompositions of Equation|6] First experimen- 
tal results are discussed in Section[3] 



3 Present Status: The Scenario Supported by Femtoscopic Data 

In the absence of space-momentum correlations, s{'p,x) would factorize into a func- 
tion of p multiplied by a function of x, and all source parameters would be inde- 
pendent of the pair momentum P and independent of particle species. Collective 
hydrodynamic-like flow is a principal source of such correlations, but not the sole 
cause. For instance, jets represent a particularly dramatic example of space-time 
correlations as the particles with high transverse momentum hadronize further away 
in the outward direction. Resonant decays also provide correlations to some de- 
gree. Given that the characteristic mean free paths at early times are much shorter 
than the size of the fireball, the language of hydrodynamics is usually invoked to 
describe space-momentum correlations for low to mid-p, particles in heavy-ion 
collisions. Generically, the femtoscopic data paint a clear scenario dominated by 
hydrodynamic-like flow. Such a source leaves several unique femtoscopic finger- 
prints. In this section, we summarize these signatures, and see how each of them 
has been observed in the data so far. 
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3.1 Longitudinal Flow: Rapidity and mj Systematics 



At ultrarelativistic energies, one expects a source extended- and with strong flow- 
along the beam axis; this is in contrast to lower energies, where three-dimensional 
radial flow is observed (l67t roSc l69h . Indeed, many models assume pure boost in- 
variance, at least at midrapidity, in which the space-time rapidity equals the mo- 
mentum rapidity. Rapidity distributions alone suggest approximate boost invari- 
ance at midrapidity, but femtoscopic probes provide stronger geometric support. 
At RHIC and SPS, a boost-invariant scenario is supported by both Yano-Koonin ve- 
locities and the Bertsch-Pratt cross terms Jt^) (c.f. Section [2.2.2| i 
within a unit of midrapi dity . At more forward rapidities, a breakdown of this sce- 



nario IS revealed (I75t 173 1 I71tl74l) below RHIC energies. 



As seen in Figure [T] the Yano-Koonin source velocity (c.f. Section l2.2.2b is al- 
most identical to that of the pion pair, for a broad range of collision energies. Fig- 
ure |2] shows the same analysis in e+ — collisions by the OPAL Collaboration at 
LEP®; here, the analysis is performed along the event thrust direction, which is 
a more natural symmetry axis than the beam direction. 



OPAL 



-, — 1 — I — 1 — I 
<kj> 

O.IOGeV E877 

O.IBGeV NA49 

0.45 GeV WA97 

0.45 GeV PHOBOS 

0.45 GeV PHOBOS 



Fig. 1 Yano-Koonin rapidity versus pair rapidity, 
for pion pairs from central collisions, measured for 
heavy ion collisions over a range of energies. Open 
(closed) symbols show results for negative (posi- 
tive) pion pairs. (Figure from JTsh ) 




Fig. 2 Yano-Koonin rapidity versus pair 
rapidity, for pion pairs from e+ — col- 
lisions at LEP, measured by the OPAL col- 
laboration. From 



Perhaps the best known, and the most misunderstood, signature of collective flow 
is known as m, scaling (F77!), which often simply refers to the fall of the Gaussian 
radius parameters with p,. However, the source of this behavior is very different 
for Rside and Rout vs. ^long- The term nit scaling refers only to the pt dependence 
of R\ong, and is inherently a relativistic effect. If one considers a purely longitudinal 
Bjorken expansion ( ItsI) where the matter thermally dissociates at a fixed proper time 
T = Vf^ —z^, particles of transverse momentum p, and zero rapidity are emitted 
from position z = Tsinhrj, where the collective rapidity is also 77 with probabiUty 
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dN 
dzd^p 



For large z, the distribution falls off exponentially with a characteristic length, 

R -^^ 



(12) 



(13) 



whereas for small z the curvature mimics a Gaussian with a characteristic radius 
parameter ^TT\) 



m,. 



(14) 



This last relation motivated the term "m, scaling", as it suggests that ^long would 
scale as 1/ y/nk. The effect is illustrated in Fig. [3] 

Indeed the pr dependence of Rio,^^ seen in Figure |6l goes as nij", where a is 
slightly larger than 0.5 as expected due to the additional effect of transverse 
flow (7^ 80; 81). Lepton collisions at LEP display an m^-dependence of the radii 
similar to that from heavy ion collisions, however; compare Figures |6] and |7] 





Fig. 3 Emission zones for higher mass (blue) and lower mass (red) particles, as seen from along- 
side the beam (left) and looking down the beam pipe (right). Due to collective expansion particles 
of a given velocity come from a subset of the overall region, referred to as the region of homo- 
geneity. Particles of higher mass, or higher transverse mass, have smaller thermal velocities and 
are thus confined to a smaller region by collective flow. For non-identical particles with velocities 
ahead of the surface, the more massive species are more confined to the surface, which leads to an 
offset 4 r of the centroid of two outgoing phase space clouds. (Figure from ffsh ) 



Thus, both the rapidity and the mj- dependence in the data support the com- 
mon assumption of boost-invariance built into several parameterized models used 
to study transverse flow, discussed next in Section [372l 
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Fig. 4 One-dimensional Gaussian radii for 
source of various particles from central Au-l-Au 
collisions at RHIC. (Figure from (T^)-) 



Fig. 5 One-dimensional Gaussian radii for 
source of various particles from e+ — col- 
lisions at LEP. From (83) 



3.2 Transverse Flow: mj Systematics and Non-identical Particle 
Correlations 

Several parameterizations have been applied to HBT data which are based on ther- 
mal emission from a radially and longitudinally expanding source dsil IssI : l84l : [ssl : 
[86l) . The term "blast wave" refers to any number of parameterizations which usually 
include a transverse radius, a maximum transverse velocity at the surface, a breakup 
time (which in a Bjorken geometry determines the longitudinal velocity gradient, 
dv/dz = 1 /ly), and a temperature and chemical potential. Parameterizations might 
or might not include a surface-diffuseness or an emission duration. Whereas most 
parameterizations assume a step-function profile, the Buda-Lund parameterization 
has a Gaussian shape and can be motivated by exact solutions to a specific solution 
for non-relativistic hydrodynamics (i83i : i84i) . Femtoscopic signatures of transverse 
flow from these parameterizations are rather generic. The consistency of the longi- 
tudinal systematics with an approximately boost-invariance scenario, discussed in 
the previous Section, motivates the assumption of such a scenario in most models. 

The p, dependence of /?side and 7?out is generated by transverse flow0, and does 
not necessarily follow any scaling formula. For instance, if one were to consider an 
infinite Hubble flow, all three radius parameters would be independent of momen- 
tum (if expressed in the pair frame) and would be equal to XyJYJm. However, if one 
cuts off the Hubble flow at a given radial size Rmax and confines emission within that 
volume, the effective source sizes are more strongly reduced for high momentum 
particles, for whom a larger fraction of the emission region would otherwise reside 



' For sources of order 1 fm, some of the pj dependence of Sout is also generated by the resonance 
halo. 
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Fig. 6 Blast-wave jSli) calculations are 
compared to single-particle spectra, ellip- 
tic flow, and Gaussian "HBT radii" from 
Au-l-Au collisions at RHIC. Data are plot- 
ted as a function of transverse momen- 
tum relative to the beam direction. (Fig- 
ure from (13)) 
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Fig. 7 Gaussian "HBT radii" from — colli- 
sions at LEP, as a function of pair transverse momen- 
tum relative to the thrust direction. From 123). 



outside /?max- The effect is indicated by the smaller homogeneity regions for faster 
(or higher mr) particles in Figure [3] 

The femtoscopic flow signature is well-described by parametric models which 
also quantitatively reproduce momentum-only observables such as pr spectra for 
identified particles and elliptic flow, as seen in Figure |6] Including realistic reso- 
nance contributions varies the parameters (e.g. temperature) of the source, but not 
the requirement of a flow-dominated scenario (1851) . The femtoscopic flow signature 
is, in fact, observed over a huge range of conditions (0). System timescales are 
observed to grow with bombarding energy (i74l) and collision centralities JHol) . In 
all cases, freeze-out temperatures of order 120 MeV and flow velocities on order of 
0.5c are indicated by the kinematic dependences of the source radii. 

A further test of the flow scenario is to look at other particle types, in addition to 
pions. Homogeneity regions have been probed for the most common particle types. 
Table [T] indicates some of the studies reported at RHIC. A flowing system would 
produce homogeneity lengths following an approximately universal dependence on 
niT, independent of particle type. As seen in Figure |4] a universal scaling is, in fact 
observed in heavy ion collisions. A similar mass dependence is shown in Figure |5] 
for lepton collisions at LEP. 

Probably the most convincing test of transverse flow is its generic signature of 
the new information accessible in non-identical particle correlations, discussed in 
Section 12] In particular, as indicated in Figure [3] not only do higher-wzj- particles 
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Table 1 A table of published or ongoing femtoscopic studies at RHIC for various particle combi- 
nations. "Traditional" identical-particle interferometry lies along the lowest diagonal line of cells. 
The table lists representative reports only; it is not exhaustive. 



come from a smaller homogeneity region than low-mj- ones, but their average emis- 
sion point lies further away in the "out" direction (ISlh . Few published data on non- 
identical particle correlations at RHIC are available, but STAR has reported the 
offset between pions and kaons (15 7k to be between 4 and 7.5 fm, consistent with 
blast wave expectations of 7 fm ( Isih . Several measurements of non-identical parti- 
cle correlations are underway at RHIC; they appear in the non-diagonal elements of 
Tabled 



3.3 Overall Source Size and Shape 

Collective and random motion (parameterized by flow and temperature, respec- 
tively) determine the femtoscopic substructure of the system created in heavy ion 
colhsions. However, they do not determine the overall space-time scale of the 
freeze-out configuration. Consider two "blast-wave" systems, with identical tem- 
peratures and flow; the second system is identical to the first, but scaled down in 
space-time. They would produce identical momentum-space observables, spectra 
and elliptic flow. Femtoscopically, the pj, y, <p and non-identical particle systemat- 
ics would be identical, except for an overall factor of two scaling in extracted length 
scales ("HBT radii" etc.). What determines the overall length scale? 

As shown in Figure |8] the CERES Collaboration studied the energy evolution 
of Rioni^R^ide fo'" Pt pious from central collisions between the heaviest nuclei; 
this quantity corresponds roughly to the overall freezeout "volume." Assuming a 



universal mean free path at freeze-out, they found ( l95h that this "volume" seemed to 
be driven solely by the final-state chemistry- i.e. the yields of pions, protons, etc- of 
the system. Above AGS energies (^/snn ^ 5 GeV), the system is pion-dominated, 
so the determining quantity is dNji/drj w dNch/drj. 
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Fig. 8 "Effective pion cross-section" (left axis) and pion freezeout volume (right axis) as a func- 
tion of collision energy, for central heavy ion collisions. From i9^ . 



Actually, the statement can be made much more generally. In Figure |9] we see 
that both Rside ™d Rumti, individually, scale with dNciJdr]. Rside is related most di- 
rectly to geometrical length, and Riong to timescale i8l|) . Rom depends on both 
space and time in a more convoluted way (9^ 81); consequently, the multiplicity 
scaling is broken for this HBT radius. Since it includes results from a range of cen- 
tralities in each collision system, Figure|9]makes clear that the chemistry-scaling is 
not restricted to central collisions. Going further, the universal -dependence men- 
tioned in Section l372l implies that the chemistry-scaling is not restricted to any given 
momentum region. Finally, inclusion of Si-induced collisions on Figure|9]shows that 
the scaling is not restricted only to collisions between the heaviest nuclei; indeed, 
preliminary data (1971 ; 1981) from the STAR Collaboration on Cu-i-Cu and Au-i-Au col- 
lisions at various energies, and even from d-i-Au and p-i-p collisions, further support 
this conclusion. 

As seen in the compilation ( l99l) of Figure [TOl a similar multiplicity dependence 
is observed in high energy hadron collisions. Since the measurements shown were 
performed with somewhat different analysis techniques and momentum acceptance, 
comparing the absolute magnitude of the radius experiment-to-experiment is non- 
trial. However, within a given experiment, the multiplicity dependence is clear 



3.4 Anisotropic freezeout configuration: systematics relative to the 
reaction plane 

Overall particle yields, rapidity distributions, pr spectra and elliptic flow (v2) are all 
measurements of the number of particles of a given momentum. So far, we have seen 
a considerable correspondence between these observables and femtoscopic ones. 
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Fig. 9 Three-dimensional "HBT 
radii" for pions from central heavy 
ion collisions over a broad range 
of energy, plotted as a function of 
charged-particle multiplicity. 



Fig. 10 "HBT radii" for pions from p — p, p — p, and a — 
a collisions at the Tevatron, SppS and ISR, as a function of 
multiplicity. From a compilation in i99l) . 



Fully integrated number yields are typically used to estimate the temperature and 
chemical potentials in a chemical equilibrium scenario. Meanwhile, the overall fem- 
toscopic scale (or volume) is driven also by the chemical composition at freezeout. 
Similarly, the shape of dN/drj reflects the degree of boost-invariance of the system, 
as do the longitudinal systematics of femtoscopic scales. Transverse momentum 
spectra carry information on explosive transverse flow, as do the pr systematics of 
femtoscopy scales. 

By far, the most intensely studied soft-sector observabl e at R HIC is elliptic 
flow- the number distribution relative to the reaction plane dlOOl) . Far more than 
the azimuthally-averaged distribution, the anisotropic part of the momentum distri- 
bution sensitively prob e s imp ortant details in transport models, such as the Equa- 



tion of State (e.g. llOlt 11021) . So far, corresponding femtoscopic studies are few. 



but they may be as important as elliptic flow measurements. Azimuthally-sensitive 
femtoscopy reveals the anisotropic shape of the emitting source at freezeout. Know- 
ing the strength of the anisotropic flow and the initial geometric anisotropy from a 
Glauber calculation, one may estimate the tir ne re quired for the system to evolve 
from its initial to final shape. Such estimates ( Il03h yield evolution timescales con- 
sistent with those estimated by R/oni; (cf. SectionITT)- 

Whereas the overall magnitude and the pj- and y-systematics of femtoscopic 
scales vary little with bombarding energy (15), the anisotropic shape varies con- 
siderably. Figure [TT] shows the world dataset of azimuthally-sensitive femtosopy 
measurements. The source anisotropy e is defined and estimated (18 ih according to 



20 



M.A. Lisa and S. Pratt 



\E89S 


CERES 


STAR 


ALICE 






▲ 








★ 






* 




99 


- — ROMD 








A hydro 


,1 , , , 


-1 





2 3 4 

1 10 10 10 10 

Vs^N (GeV) 

Fig. 11 The transverse spatial freezeout anisotropy e as a function of collision energy, for mid- 
central (~ 10 — 30%) heavy ion collisions. Round sources correspond toe = 0; e > indic ates 
an out-of-plane-extended source. Measurements at the AGS jl04) . SPS jlOSi) and RHIC il06h are 
comp ared with RQMD transport model calculations at the low energies; hydrodynamic calcula- 
tions ilOTh at RHIC and the LHC are indicated. 



Ky — Kj. 

Rl + Rl 



'R 



(15) 



.v,0 



where Rx (Ry) represent source dimensions in (out of) the reaction plane. R^.q and 
R^.2 ^re the azimuthally average of the Rside "HBT radius" (c.f. Section |2. 2. 2l i and 
the strength of its oscillation with azimuthal angle. 

At t he A GS, the emission zone is almost as anisotropic as the nuclea r ove rlap 
region (11041) . not surprising since elliptic flow vanishes at these energies (llOSh . As 
both the evolution duration and strength of anisotropic flow increase with ^/snn, the 
source becomes more round, and is even expected to become in-plane-extended at 
the LHC, according to hydrodynamical calculations ( 107 : 86l) . Taken at face value, 
the data suggest a non-monotonic energy dependence, perh aps w ith a dip at SPS 
energies, where there are other hints of threshold behavior (11091) . However, much 
more systematic data is needed to confirm this intriguing energy dependence. 



3.5 Resonance decay contribution: non-Gaussian "imaging" fits 

While Gaussian radii from fits to correlation functions capture the bulk spatial 
scales, the emitting system clearly has a richer structure. In particular, one expects 



an exponential component to the source distribution (a "halo" (l58l) ) due to long 



and medium-lived resonance decays. As discussed in Section 12.2.31 in "imaging" 



fits (l65t l63l). a spline-based separation distribution S{r) is varied such that the cal- 
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culated moments of the correlation function (c.f. Equation |6l) match the measured 
ones. 

Separation distributions so extracted are shown in Fi gure[T2l from two-pion cor- 
relations from central Pb+Pb collisions at the SPS Exponential-like tails arise 
in all three dimensions, and are due mainly to resonant decays, but along the lon- 
gitudinal direction are also influenced by the exponential tail coming from boost- 
invariant dynamics as illustrated by Equation [12] The Therminator variant of the 
blast wave model also incorporates resonant decays (fssl 'si'), and represents the 
state of the art in such models. On the left side of Figure [12] is shown Therminator 
calculations of the Gaussian "HBT radii," compared to STAR measurements. 

Many fewer resonances are created at the much lower energies (^/snn ^ 3 GeV) 
used in AGS experiments. Thus, it is not surprising that the long-range exponential 
tail of S{r) is not observed there (l66l) . 



3.6 Putting it all together: factorization and the overall picture 

The decades-long femtoscopic program in heavy ion collisions has systematically 
turned most of the available "knobs" to map out the landscape. These include: 
initial-state conditions, such as ion size and collision centrality and energy; kine- 
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matic knobs such as pair rapidity and azimuthal angle; and femtoscopy between dif- 
ferent particle types. Overall, the data show a remarkable "factorization." Denoting 
a generic femtoscopic observable (e.g. Rumg) by 3, all data indicate, approximately 

J(\/iAW,A,B, |b|,^,y,OT7',mi,OT2) « (16) 
Rg {dNn/dri,dNp/dri,. . .) ■ F„ . . .) ■ Fk {y,mT ,mi,m2) , 

where A and B represent the size of the colliding nuclei, mi and m2 are the masses 
of the two particles being correlated, and dNi/drj is the rapidity density of particle 
i. 

According to "Equation" [16] the overall, dimensionful scale of the system is 
determined by the collision conditions ^/snn. A, B and |b|. Going further, there is a 
"universality" in this dependence: it is not these initial-state conditions per se, but 
the final chemistry of the collision (dNi/drj) which determines the overall factor 7?^. 

The transverse and longitudinal collective dynamics of the system is character- 
ized by the dimensionless function F^, which essentially does not depend on colli- 
sion conditions or chemistry. Blast-wave type collectivity leaves unmistakable sig- 
natures on this function, which is supported by all observations. Indeed, as suggested 
by Figure]?] the dependence of F/^ on the masses of the correlated particles only ap- 
pears when mi ^ m2, for which a shift in the average emission point may occur; c.f. 
end of Section [J!2] 

The azimuthal dependence of the source geometry is singled out, since its depen- 
dencies have not been fully mapped. Clearly, there is at least some depe ndence on 



both in itial-state and kinematic variables: ^/snn (e-g- FigurefTTT). pr (Il06l) . an d cen 



trality ( 1061 ; 105 ). Given its clear connection to elliptic flow, evolution time ( 103 ) 



significant energy evolution, and predicted sign change at the LHC, further femto- 
scopic analyses relative to the reaction plane are of high priority. 

When 2 = Rout, the factorization of Equation[T6]is somewhat violated, as shown 
in Figure |9] As with the azimuthal dependence of femtoscopy. Rout clearly mixes 
both space and time ( l96l) , which scale somewhat differently with initial conditions. 

All available femtoscopic data point to a system dominated by collective motion: 
approximate longitudinal boost-invariance, strong transverse flow of order 0.5c, a 
random kinetic component characterized by a "temperature" ~ 0. 1 — 0. 15 GeV. The 
dimensional scale of the system at freezeout is determined by a universal hadronic 



mean free path ~ 1 fm (I95h ; its transverse size is typically 2x the initial overlap 
zone (l60l) . though feed-down from resonance decay add a long-range exponential 
tail. 

The timescales of the collision are short. The "HBT radius" Ri„ng and the differ- 



ence between the initial and final shape of the source are both consistent (11031) with 
an evolution time^ 10 fm/c for central heavy ion collisions. The duration of particle 
emission is short- on order 3-4 fm/c for central collisions. Again, resonance decays 
add a tail to the time distributions. 

These conclusions are based largely on blast-wave fits to available data. Although 
such models have numerous parameters, three things should be kept in mind. Firstly, 
the parameters do represent a minimal set (c.f. e.g. fs]]). Spectra as a function of 
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mass and /jj are determined by one parameter each for flow and temperature. Once 
these parameters are set, they leave unavoidable fingerprints on the function 7?^ from 
"Equation" [161 they completely determine the momentum and mass dependence of 
identical and non-identical correlations alike. These fingerprints cannot be changed 
by the remaining parameters, which specify the size and timescales of the system 
and are reflected in the data by the function Rg. 

Despite the remarkable success of blast-wave models, one wonders whether other 
factors might imitate flow. For heavy ion collisions, an obvious candidate is cool- 
ing. If emission occurs over a more extended time, but with a falling temperature, 
one would also expect radius parameters to fall with pt. However, cooling can only 
be important if the duration of the emission is significant, and extended durations 
are incompatible with the experimental observation that /?out/^side ^ 1- Thus, it be- 
comes difficult to offer any plausible counter-scenario, to compete with the picture 
of a collision dominated by a rapidly developed collective radial flow, followed by 
a sudden dissolution. 



3. 7 What about small systems? 

The study of ultrarelativistic heavy ion- as opposed to hadron or lepton- collisions 
is driven by the desire to create a bulk, self-interacting system whose properties 
and response may eventually be linked to an equation of state (EoS). The EoS is 
not characteristic of collisions, per se, but is a more fundamental object, directly 
reflecting the symmetries and properties of non-perturbative QCD. The use of heavy 
ions is based on their size- a "large" system (whose scale is much more than the 
mean free paths of particulate degrees of freedom) is more likely to equilibrate. 

The generation of a deconfined state of matter- in which colored degrees of free- 
dom are relevant over length scales greater than a hadron radius- would also seem 
to call for "large" colliding systems. Apparently, particle collisions need not apply. 

As we have discussed, femtoscopic measurements in heavy ion collisions all 
point to a flow-dominated scenario consistent with all other soft-sector observables. 
This flow has become both the greatest evidence of the creation of a bulk system at 
RHIC and a tool to probe its nature. Its signature (inscribed in the function Ri^) is 
seen for a huge range of initial conditions. 

Indeed, as we have seen in Figures |2] |5] |7] and [TOl femtoscopic systematics- 
which seem to make much sense in heavy ion collisions- are found also in hadron- 
hadron and e+ — collisions (c.f. |99l for a compilation). Due to technical differ- 
ences in femtoscopic studies in particle and heavy ion measurements, direct compar- 
isons between these measurements is difficult. However, the first apples-to-apples 



comparison of pion femtoscopy from p + p and Au +Au collisions at RHIC (l97h . 
using identical methods and detection conditions, reports an essentially identical 
pr-dependence of R^, and clear multiplicity dependence of Rg, consistent with pre 
vious measurements i9. 
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Femtoscopic measurements are not the only soft-sector observables displaying 
very similar systematics in heavy ion and particle collisions. In hadron collisions 
at the highest energies, the mj d i stributions for heavier particles are less steep 



that those for light particles (111 It Ill2t Ill3h . just as they are in heavy ion col- 
lisions (c.f. Figure |6]l. Recent studies ( 114l) suggest that the observed breakdown 
of this similarity at low energies (11 15h may be due to simple energy and mo- 
mentum conservation effects. Similarly, conservation of discrete quantum numbers 
(strangeness, baryon number) leads to differences in particle yield ratios measured 
in heavy ion an d — collisions; but when these are accounted for, the trends 
are very similar 161) . In fact, thermochemical fits to yields from p + p collisions at 
= 200 GeV (11 17b are essen tially identical to similar fits to AuH-Au collisions at 
the same energy (e.g. Fig 12 of ll 181) . 

The data obviously beg the question: if the soft sector in one collision system 
appears identical to that in another, is there any firm justification to assume that the 
physics driving these observables is different? 

In considering this question, some important points may be kept in mind. One 
may imagine that string-breaking or mini-jet dynamics dominate the correlations 
in hadronic collisions, so that perhaps the same mechanisms drive 7?^ in heavy ion 
collisions. However, the combined observations of Sections |3.2| and 13.31 directlv 
contradict this. A minijet or string has a characteristic scale, presumably ^ 1 fm. 
The growth of the homogeneity region with multiplicity (e.g. Figure |9]l might be 
understood in terms of correlations between particles from different minijets. How- 
ever in this case, there would be very little p 7- -dependence of 7?^; particles with high 
momentum would arise from all regions of the "whole source," as would slower par- 
ticles (13 8h . Similarly, the suggestion (l82h that the Heisenberg uncertainty principle 
explains the space-momentum correlation in small systems cannot apply equally to 
heavy ion collisions; there is no good reason to imagine a "scaled uncertainty prin- 
ciple" in which pr ■ R fa xK where x grows to ~ 6 for the most central heavy ion 
collisions. 

The striking resemblance of hadronic and nuclear collision systems deserves 
much further theoretical focus and experimental exploration at RHIC and the LHC, 
including measurement of non-identical particle correlations in small systems. 



4 Using Transport Models and Femtoscopy to Determine Bulk 
Properties 

Any reasonable inference of bulk properties of matter such as the equation of state 
or viscosity from data requires detailed analysis with sophisticated models of the 
entire evolution. Thus one can not understand static properties of bulk matter with- 
out simultaneously understanding the dynamics. The closest thing to an exception 
to this rule concerns the entropy. Using a combination of femtoscop y and spectra, 
one can measure the net entropy per unit rapidity of the final state dl 19t) . This is 
accomplished by first extracting the phase space density by dividing the spectra by 
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the femtoscopic volume. If one assumes that entropy is conserved, one could in- 
fer a point on the equation of state if one knew the energy density and volume at 
some early time where the matter has thermalized but before transverse flow devel- 
oped. At a time of T = 1 fm/c, one can state the volume of the source with some 
confidence, Y = ttzR", where R is the initial radius of the fireball. However, the 
net energy is less certain. The transverse energy per rapidity can be measured for 
the final state, but neglects the energy associated with longitudinal thermal motion, 
and also neglects the longtudinal work done by the expansion between T = 1 fm/c 
and breakup. Nonetheless, the final state entropy, ^ 4450 ± 10% units for central 
130A GeV AuH-Au collisions, was in the range expected from lattice gauge theory. 
Of course, this inference of the initial entropy is modified if there is large entropy 
generation or if the matter has not thermalized at T = 1 fm/c. 

When the first pion correlation measurements at the SPS were performed by 
NA35 in 1988, it was hoped that femtoscopy would reveal long lived sources, 
proving the existence of a strong first-order phase transition. The reasoning behind 
the expectations was based on model calculations assuming matter had completely 
stopped, and thus would not reaccelerate except for the pressure. For matter at the 
border of the mixed and QGP phases, the ratio of pressure to energy density, P/e, 
would be a minimum, and one would observe a maximum of the lifetime. For simple 
bag model equations of state, this ratio was expected to fall to approximately 3/40 
at energy densities near 2 GeV/fm^. If the matter had indeed stopped, and if the bag 
model equation of state were true, the system would require more than 50 fm/c to 
dissolve Ins S). 

Such dramatic lifetimes were not observed at the AGS and SPS, even though the 
expected energy density was attained somewhere in the lower SPS region. Two main 
reasons contributed to the failure. First, as expected, the stopping stage was trans- 
parent which mean the system would dissolve due to longitudinal expansion which 
did not rely on pressure or transverse flow. Nonetheless, one would have expected 
lifetimes near 20 fm/c even with the transparency. The fact that it appears that disso- 
lution times were closer to 10 fm/c is mostly due to the fact that the true equation of 
state significantly varies from that of the simple bag model equation of state. Cur- 
rent lattice calculations show no first-order phase transition and show P/e staying 
more than a factor of two above the bag-model minimum (11211) . Additionally, the 
critical energy density is reached for the lower SPS range, where the non-negligible 
net baryon density which might provide extra repulsion. 

Energy densities at RHIC surpass those of the SPS, and at 1 fm/c are in the 
neighborhood of 10 Ge'V/fnv', roughly a factor of 5-10 higher than those where one 
might expect a minimum in P/e. As one would expect, femtoscopic inferences of 
the dynamics at RHIC show somewhat more explosive collisions than at the SPS or 
AGS. Although the trends with beam energy were qualitatively in line with expec- 
tations, quantitative predictions with the hydrodynamic models of the early RHIC 
era ( Il23l : ll02h . which had been so successful in describing elliptic flo w ob servables, 
provide poor reproductions of the data. The failure of the model in ( Il23h was par- 
ticularly disquieting given that it incorporated a cascade afterburner which should 
provide an accurate description of the breakup stage. In these models, the character- 
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Fig. 13 Gaussian radii reflecting spatial sizes of outgoing phase space distributions in three direc- 
tions: Rout, Rside and Sioiia. Data from the STAR collaboration (red stars) are poorly fit by a model 
with a first-order phase transition, no pre-thermal flow, and no viscosity (solid black squares). 
Correcting for all those deficiencies, and using a more appropriate treatment of the relative wave 
function in the Koonin equation JTslilsl) brings calculations close to the data (filled black circles). 
The sequential effects of including pre-thermal acceleration (open blue squares), using a more real- 
istic equation of state (open green diamonds), and adding viscosity (open cyan triangles) all make 
substantial improvements to fitting the data. An improved relative wav e fun ction yielded modest 
improvements (compare open cyan triangles to filled black circles). See jl22l) 



istic time for emission (related to /?iong) of particles seemed to be « 15 — 20 fm/c, 
rather than the < 10 fm/c one would infer from blast wave analyses. Furthermore, 
the duration of emission was significantly longer in the models than in the blast- 
wave analyses of the data, as the models predicted ^om/^side ^ 1-5, as opposed to 
the experimental value of « 1 . In contrast, cascade models, which many would have 
expected to fail given the lack of important f eatures in the eq uation of state, provided 
a significantly better description of the data (1 124c 1 1251 : 1 1261) . The effective equations 
of state for the cascade models was stiffer than those used in the hydrodynamic mod- 
els, with some softness provided by resonances. A cascade model written without 
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the inclusion of resonances (Il27l) significantly under-predicted experimental source 
sizes. 

The disparate array of models mentioned above appeared to provide a consis- 
tent trend, in that models with stiffer equations of state provided more explosive 
dynamics, which lead to smaller Gaussian source parameters and smaller ^out/^side 
ratios. However, one must remain wary of several significant shortcomings of the 
models, which might significantly affect any conclusions generated by comparison 
to femtoscopic data. The following list of model features and parameters have all 
been shown to significantly influence femtoscopic results. 

• Stiffer equation of state: A stiffer equation of state signifies a reduced entropy, 
which for a fixed spectra gives smaller femtoscopic volumes. Additionally, a 
higher pressure leads to more rapid expansions which gives smaller /?out/^side ra- 
tios. Whereas calculations from earlier in the RHIC era used first-order equations 
of sta te ( Il23h . incorporating a finite speed of sound, as seen in lattice calculations 
( 121 ), in the region near Tc can lower the /?out/^side ratio at the 10% level (1128|). 

• Improved modeling of non-equilibrium hadronic chemistry: Significant acceler- 
ation occurs during the hadronic phase, during which chemical abundances tend 
to be over-populated. This lowers the pressure at the periphery which helps a 
pulse form which leads to a mo re sim ultaneous disassociation and an /?out/^side 
ratio which is lower by ~ 10% (Il29t) . 

• Longitudinal Acceleration: Nearly all hydrodynamic codes used for femtoscopic 
comparison have used boost invariance, where all fluid elements have zero ac- 
celeration in the longitudinal direction, to simplify calculations. Including longi- 
tudinal acceleration breaks the simple relation between the longitudinal velocity 
gradient and the lifetime. This shortcoming might lead to models overstimating 
^long, and for parametric blast-w ave analy ses could lead to an underestimate of 



true colUsion Ufetimes by 5-10% (ll30l:ll3lh . 



Shear and bulk viscosity: At early times, where the velocity gradient is large and 
longitudinal, shear alters the isotropy of the stress-energy tens or by lowering the 
effective longitudinal and increasing the transverse pressure (Il32h . This accel- 
erates the transve rse acceleration and can lead to a lowering of the Rowi/RsHe 
ratios by 10-15% ( Il33t 11341 11281) . Bulk visco sity leads to a lowering of the ef- 



fective pressure for energy densities near Tc ( Il35t 11361) which strengthens the 



pulse. However, for the periphery of the collision, the pressure is unchanged and 
the pulse ten ds to largely dissipate by breakup. The effect on source radii is then 
rather small (Il28h . 

Initial energy density profile: The profile is affected by the scheme for screening 
collisions. The three most often-used schemes are the w ounded-nucleon model 
( 137 : 138 ). color-glass motivated prescriptions (139; 140l) and collisional scaling. 
Of these three prescriptions, collisional scaling gives the most compact fireball 
and the wounded-nucleon model leads to the least compact. More compact fire- 
balls are m ore explosive and result in a lowering of /?out/^side ratios by 5-10% 

(EH nil. 

Pre-equilibrium transverse acceleration: At early times, T < 1 fm/c, the departure 
from equiUbrium are likely so large that modifications to the stress-energy ten- 
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sor can not be described by Navier-Stokes concepts. If the matter is dominated 
by weakly longitudinal color fields the trans verse pres sure is significantly higher 
than that of an equilibrated relat ivistic gas (Il42tll32h . This can again lower the 
^out/^side ratio by ^ 10% (Il28b . Similar conclusions have been reached by in- 
serting an initial velocity profile into ideal hvdrodvnamic s (1 1431) and adding a 
strongly repulsive potential into a hadronic cascade model (11441) . 

None of the items in the above list was able to account for the initial ^ 50% dis- 
crepancy in the /?out/^side ratio, by themselves. But, by incorporating all the above 
features into a model, the source radii were reasonably well reproduced ( 11281) . Fig- 
ure[T3]compares experimental source radii to those extracted from an ideal hydrody- 
namic model with a first-order phase transition and no pre-equilibrium acceleration. 
Successively correcting each of the se de ficiencies, one finds increasing agreement 
with the data. For more details, see (Il22h . 



Within the last year a consensus appears to be developing (e.g. Il41h that a sud- 
den impulse to the flow at early time is required to explain Gaussian source-size 
parameters. However, the relative contribution from shear and the non-equilibrium 
phase is not well known. Additionally, it is difficult to disentangle the contribution 
from adjusting the equation of state. It is hoped that a more global analysis, includ- 
ing spectra, correlations and flow observables, will be able to better distinguish and 
determine individual aspects of the physics. 



5 Femtoscopic Expectations at the LHC 

For soft sector observables in heavy ion physics, long-term baselines have been es- 
tablished over a large energy range. Prior to first data at RHIC, it was commonly 
speculated (and hoped) that large deviations from these systematics (e.g. n^/n"^ 
ratios, sidewards flow, strangeness enhancement, total multiplicity) would signal 
clearly the qualitatively different nature of the system created the re (1261) . In fem- 



toscopic systems, rather generic arguments led to expectations (I25t Il45h of a rapid 
increase, with y^s^N, in the pion "HBT radii" Roui and Rio„i;, reflecting relatively 
long timescales of the transition from deconfined QGP to confined hadronic matter. 

Such dramatic speculations are largely absent today in anticipation of LHC colli- 
sions. Soft-sector, global observables at RHIC are only quantitatively different than 
they are at lower energies. Even in the high-p^ sector, where jet suppression and 
partonic energy loss measurements have generated huge excitement, energy scans 
at RHIC reveal that the data indicate more of an evolution than a revolution. 

The need for such systematic comparisons is not unique to RHIC, but has been 



a generic feature of all heavy ion programs (1231: 1241) . from low-energy facilities like 
the NSCL (Michigan State), to progressively higher-energy facilities at SIS (GSI), 
file Bevatron/Bevalac (Berkeley Lab), AGS (Brookhaven), and SPS (CERN). The 
nature of heavy ion physics is such that little is learned through study of a single 
system. Discoveries via sharp jumps a la superconductivity are not our lot. The real 
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science behind heavy ion measurements (at very high energies as well as at much 
lower ones) lies in understanding the details of the data. 

With this in mind, what may await us in heavy ion collisions at the LHC? After 
discussing a simple multiplicity-based extrapolation of existing systematics, moti- 
vated by the factorization of Section [331 we summarize predictions of Boltzmann 
and hydrodynamical transport calculations. We will find that, in both cases, the 
much greater explosive flow predicted at the LHC may lead to new effects which 
break the simple extrapolation. 



5.1 Nothing New Under the Sun (NNUS) Scenario 



As discussed in Section I3T61 most femtoscopic systematics so far measured may be 
summarized in terms of a factorization. While the kinematic dependences may be 
described by a single dimensionless function Z^., there is an overall scale Rg which 
grows roughly as the cube root of the charged particle multiplicity. 

While the factorization is somewhat broken when 2 = and perhaps by 
azimuthally-sensitve femtoscopy, the overall multiplicity is probably our best, zero- 
new -physics guide to simple extrapolation of femtoscopic trends measured over two 
orders of magnitude in ^^s^n and from from the lightest {p + p) to the heaviest 

{Pb + Pb) systems. Figure|9]suggests a simple fovm Rg{dN /dy) °^ {dN /dy)^^^; this 
ignores the finite offset ~ 1 fm when extrapolating dN/dy 0, but this is negligi- 
ble for high multiplicity. This relation may reflect that a constant freezeout density 



drives the femtoscopic scales (l95h . though this neglects any dynamic effects. As- 
suming that this simple proportionality continues, then, determining femtoscopic 
expectations boils down to anticipating the multiplicity at the LHC. 

A naive extrapolation (il46i : il47i) of systematics suggests that dN/dy at the LHC 
will be 60% larger than that observed at RHIC. Thus, the zeroth-order expectation is 
that length scales at the LHC will be 17% (L6^''^ = L17) larger than those measured 
at RHIC, for all kinematic selections and particle species, according to EquationfTSl 

Going beyond si mple extrapolation to include a physical picture, saturation- 



based calculations (11481) give much higher multiplicity- roughly triple that at 
RHIC. This leads to expectations of length scales 45% higher than those at RHIC. 
Thus, Riong for pions at midrapidity and low pj in central colUsions would be 
1.45x7fm==10fm. 

Multiplicity predictions based on Boltzmann/cascade calculations can be sig- 
nificantly higher yet. Selecting two for which femtoscopic pred ictions also exist 



(Section|5]21i, A Multi-Phase Tran sport (AMPT) calculation (11491) and the Hadronic 
Rescattering Model (HRM) dlSOl) predict 5 x and 7 x RHIC multiplicity, respec- 
tively. Thus, femtoscopic scales at LHC may be as much as 90% higher than at 
RHIC. Depending on the final-state interaction which produces the two-particle cor- 
relation function, measuring length scales of ^ 15 fm may challenge experimental 
two-track resolutions. For t wo-pi on correlations, such scales are within the capabil- 



ities of the ALICE detector (Il51b . 
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Pion freezeout time distribution 
••^ rescattering, b=8 fm centrality 



LHC Pb+Pb 
RHIC Au+Au 



rcscailcring, b=8 fm centrality 



. ' * 



Fig. 14 The freezeout time distribution from Fig. 15 HBT radii from fits to pion correlation 
the H adron ic Rescattering Model of Hu- functions fro m the Hadronic Rescattering Model 
manic flsoh for RHIC and LHC conditions. of Humanic flsoh for RHIC and LHC condi- 
tions. 



5.2 Boltzmann Transport Calculations 

More interesting than simple scaling relations are models with real physics and dy- 
namics, such as transport calculations. Boltzmann/cascade transport models gener- 
ally reproduce "HBT radii" at RHIC better than do hydrodynamic calculations ( isj). 
The reasons behind this include different physics in the models, a more detailed 
description of the kinetic freezeout, and the use of more appropriate methods of 



calculating the radii (Il52h . Predictions of pion HBT radii with each of the transport 
calculations discussed in Section lSTI reveal predictions more subtle than the simple 
multiplicity-scaling discussed above. 

For an infinite and boost-invariant system (only an approximation of reality, of 
course), the longitudinal HBT radius Riong is proportional to the system evolution 
time (i.e. between interpenetration of the ions and kinematic freezeout of the prod- 
ucts) (3^ 15). Naturally, this is not a unique, system-wide time, but a distribution. 



An example is seen in Figure [141 in which the pion freeze-out time distribution 
for collisions at RHIC and LHC are compared in the HRM calculation. The LHC 
timescales are roughly double those at RHIC. Although HRM is not explicitly a 
boost-invariant model, we see in Figure[T5]that Rumg reflects this timescale increase, 
roughly doubling when the energy is increased from RHIC to LHC energies. The 
~ 70% increase in Riont: is roughly consistent, then, with expectations from both 
a timescale and from the multiplicity-scaling point of view. This is certainly not a 
coincidence; at the LHC, the particles simply require more time to get isolated from 
their more numerous neighbors. 

On the other hand, there is more going on. The predicted /77-dependence of both 
Riong and Rside are steeper at the LHC than at RHIC. Also, the increase in Rsuie is 
significantly less than 90%. Both of these effects are consistent with a freezeout 



scenario with significantly increased transverse flow (1811) . Indeed, transverse mo- 
mentum distributions predicted by HRM are significantly harder (less steep) at the 
LHC than those at RHIC. Since the pr dependence of HBT radii (iTsl) and spectra 
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in the soft sector are observed to change very Httle between ^/snn = 20 200 GeV, 
it will be interesting to see whether this trend is broken at the LHC, as predicted by 
HRM. 

The HRM model is a deliberate effort to use the simplest (often criticized as sim- 
plistic) physics picture, free of novel phases like QGP. It is a pure hadron-based 
transport calculation, thou gh th e initial conditions may be taken from Pythia or 
Saturation-based scenarios (11501) . On the other side of the "simplicity spectrum" is 
AMPT, an attempt to describe the various st ages of the system's evolution in terms 
of the most appropriate model for that stage ( Il24l) . 

Similar to HRM, AMPT predicts stronger transverse flow at the LHC, as com- 
pared to RHIC, leading to steeper pr-dependence of HBT radii. In terms of scale, 
the transverse (longitudinal) radii are predicted to increase by 10% (30%). This is 
more modest than the predictions of HRM (30% and 70%, respectively), and much 
more modest than pure-multiplicity scalings of Section ISTI 

Thus the dynamical physics, in these models, lead to expected details signifi- 
cantly beyond simple extrapolation of lower-energy results. 



5.3 Hydrodynamical Calculations 

As mentioned, hydrodynamical models tend to reproduce femtoscopic measure- 
ments more poorly than do Boltzmann/cascade calculations. On the other hand, they 
have enjoyed huge success in reproducing momentum-space observables such as el- 
liptic flow. Furthermore, the conditions at LHC are likely to provide an even better 
approximation than at RHIC to the zero-mean-free-path assumptions of ideal hydro- 
dynamics. Finally, the direct connection between hydrodynamics and the Equation 
of State of strongly-interacting matter (color-confined or not) remains a compelling 
reason to explore soft-sector, bulk consequences of the model. 



5.3.1 Source Length Scales 

Recently, Eskola and collaborators (Il53b coupled a pQCD-nsaturation-based predic- 
tion for initial conditions at LHC to their 1 + 1 -dimensional hydro calculation. The 
Equation-of-State featured a first-order phase transition between an ideal QGP at 
high temperature and a hadron resonance gas at low temperature. 

As shown in Figure [16] the initial energy density at which hydrodynamics is as- 
sumed to take over expected to roughly an order of magnitude larger at the LHC 
than at RHIC, due both to increased gluon production and to shorter system forma- 
tion (thermalization) time To at the higher energy. Since the initial transverse scale 
changes only little, the pressure gradients will likewise be higher at LHC, lead- 
ing to increased transverse flow. These effects place competing pressures on the 
space-time evolution of the system, and on the femtoscopic scales at freezeout, as 
discussed below. 

The increased energy density (directly associated with entropy density and thus 
multiplicity) tends to produce longer timescales at the LHC. Longitudinal expan- 
sion tends to cool the system towards freezeout conditions. However, especially at 
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Fig. 16 Initial energy den- Fig. 17 Left: freeze-out space-time hypersurface in colli- 
sity distribution in the transverse sions with finite impact parameter calculated by Heinz and 
pl ane, calculated by Eskola et Kolb |lO"^ ). The "IPES" calculation is an estimate of the sys- 
al ll53h . tem created at the LHC. The impact parameter is defined to 

lie in the x-direction. Right: time-,z- and momentum-integrated 
freeze-out shapes in the transverse plane. 



the LHC, the large transverse flow generat ed by the intense pressure gradients can- 
not be ignored. Eskola and collaborators (Il53h . estimate that the time required to 
cool from the maximum energy density (at r = in Figure[T6lto the critical energy 
density (£,. = \ .92,GeV / fm^) would be 6 fm/c (20 fm/c) at RHIC (LHC), due to 
longitudinal expansion alone. However, when transverse dynamics are included, the 
cooling times become 5 fm/c (7.5 fm/c) at RHIC (LHC). The evolution time to kine- 
matic freezeout- say until T w 140 MeV- is Tq ~ 12 — 14 fm/c in both cases; this 
is the timescale most directly probed by femtoscopy. This is dramatic- the effect of 
transverse flow on cooling timescales can almost be neglected at RHIC, while it is 
dominant at the LHC. This is reminiscent of the cascade calculations discussed in 
Section|5]2) the much stronger flow may well lead to deviations from the trends (e.g. 
P7--dependence of pion HBT radii being independent of ^/Jnn) established so far at 
lower energy. This aspect of the NNUS scenario may finally be violated. The qual- 
itative difference is apparent from the freezeout hypersurfaces at RHIC and L HC, 
shown in Figure [17] The Figure is from a calculation by Kolb and Heinz (Il07h . but 
is similar to Eskola's. 

It would be very interesting to know whether the other aspect of NNUS, namely 
the multiplicity scaling shown in Figure|9] is satisfied by the hydro models. Unfortu- 
nately, Eskola did not calculate pion "HBT radii," and Heinz and Kolb did not calcu- 
late multiplicity, so a consistent estimate of the scaling cannot be checked. However, 
the former predict that the multiplicity at LHC will be approximately triple that at 
RHIC, corresponding to a 40% increase in HBT r adius under NNUS scaling. The 
transverse radii predicted by Heinz and Kolb (1 1071) incr ease by a similar amou nt. 

Recent hydro calculations by Sinyukov et al (11541) and Chojnacki et al (Il55l) 
report predictions similar to those of Boltzmann calculations from the previous Sec- 
tion, again due to the much larger flow. In particular, single-particle spectra are 
predicted to be much shallower at LHC than at RHIC. Rumi; is predicted to grow by 
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^1 ong 


NNUS 


17% 


17% 


17% 


HRM tl5Q) 


45% 


45% 


60% 


AMPT 11491) 


10% 


10% 


30% 


Eskola liydro (153") 


40%? 


40%? 


40%? 


Heinz hydro ( 156) 


30% 


30% 


0% 


Sinyukov hydro (154^ 


15% 


25% 


25% 


Chojnacki hydro <155|) 


25% 


45% 


25% 



Table 2 Rough indications of the increase in pion HBT radii at the LHC, relative to RHIC, pre- 
dicted by transport models. Since the models predict a stronger kj dependence at the LHC, these 
numbers will vary with kj. The numbers listed are taken at low kj (~ 150 MeV/c). Our estimate 
for Eskola's hydro is based only on simplistic multiplicity scaling; actual HBT radii have not been 
calculated in this model. 



about ^ 70%. Rout (Rside) should increase ~ 15% (~ 40%), and fall more steeply 
with pt at the LHC. These seem to be rather generic predictions from the transport 
models. They represent a deviation from the almost universal 797- -dependence of Rj^ 
(c.f. Section [l!6] l measured thus far, and thus will be an important test of a much 
more explosive source predicted at the LHC. 

While the hydro-based predictions Frodermann et al (fl56|) agree with most of 
those discussed above, they are unique in that no change at all is expected for Rumt;- 

In Table |2] we summarize the approximate increases expected for the HBT radii 
at the LHC, relative to those at RHIC. One must keep in mind, however, that the 
models predict a steeper fcj- -dependence of the radii at the LHC, so these numbers 
are taken only at low pT- 



5.3.2 Source Shape 

Azimuthally-sensitive pion interferometry- the measurement of spatial scales as a 
function of emission angle relative to the reacti on pl ane- probes the shape of the 



freezeout configuration, in addition to its scale (11571 : |48t 1811) . At finite impact pa- 
rameter, both the spatial configuration of the entrance channel and the resulting 
momentum distribution in the exit channel are anisotropic. In particular, the initial 
state is spatially extended out of the reaction plane, and the resulting flow is stronger 
in the reaction plane (elliptic flow V2 > 0). 

Due to the preferential in-plane expansion, as the system evolves the spatial 
configuration should become increasingly in-plane extended (equivalently, decreas- 
ingly out-of-plane extended). Thus, knowledge of the entrance-channel shape (e.g. 
though Glauber model calculations) and measurement of the exit-channel shape 
(through femtoscopy) provide "boundary conditions" on the dynamical spacetime 
evolution of the anisotropic system, and probe the evolution timescale. The ex- 
tracted timescale is model-dependent, requi ring in principle a detailed time evolu- 



tion of the flow. However, a simple estimate (Il03h of the timescale extracted through 
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601) to azimuthally- 



shape measurements and that extracted from blast-wave fits (jSl 

integrated HBT radii are roughly consistent. 

Measurements of the freezeout shape at the AGS (Il04 . SPS ( 1 105b and RHIC (fl58 : 

Il06l) indicate an out-of-plane-extended configuration. Consistent with the fact that 
preferential in-plane expansion (i.e. elliptic flow) is stronger at RHIC, the config- 
uration at the higher bombarding energies is rounder than at the AGS as shown in 
Figure nn 

The anisotropic shapes and corresponding azimuthally-selected HBT radii have 
been calculated in Boltzmann and hydrodynamical models. At the AGS, the trans- 
port code RQMD (159) reproduces the overall scale ( l42h and the anisotropy ( 148 



104 ) of the source reasonably well, as shown in Figure[TT] At RHIC, two groups ( 1107 
155h have reported that 2 + 1 hydrodynamical code s reproduce the shape quite well. 



At the LHC, the hydro calculations predict (llOTt l86l) a qualitative change in the 
freezeout distribution. As shown in the right panel of Figure [TT] the source is ex- 
pected to evolve to an ;«-plane configuration (e < 0). This, and the increased overall 
size, arise due to the much larger flow at the LHC, together with the increased evo- 
lution time. Geometrical and dynamical effects combine non-trivially to generate 
the radius oscillations, but at very low pr geometry dominates, and one may in fact 
measure an inversion of the sign of the oscillations ( 107 : 156h . To s ee this dramatic 
effect in ALICE, pion tracking with the Inner Tracking System (Il5lh will be crucial. 



5.4 Proton Collisions 

Before heavy ions are accelerated at the LHC, proton collisions at y/s = 1 .4 TeV 
will be measured. While the thrust of the p + p program is towards Higgs physics, 
it is well-recognized that p + p collisions serve as a valuable reference to heavy 
ion analyses in the "hard" (high-pj-) sector, where one looks for the effects of the 
medium on particles coming from well-calibrated fundamental processes. 

Soft-sector analyses, too, should be performed for systems from the smallest to 
the largest, and the results compared. Since such analyses are assumed to measure 
the bulk properties, one might well hope for qualitative differences when comparing 
results for p + p to Pb + Pb collisions. 

While pion HBT measurements have been common in both the high-energy and 
heavy-ion communities for many years, a direct "apples-to-apples" comparison be- 
tween results from A+A and p + p collisions has not been possible until very re- 
cently. The STAR Collaboration at RHIC has reported the first direct comparison 
of pion HBT radii in AuH-Au, CuH-Cu, dn-Au, and pH-p collisions, using the same 
detector, same energy, identical techniques (event mixing, etc) to create the cor- 



relation function, identical coordinate systems and identical fitting techniques ( 1971) . 
Remarkably, Gaussian fits to the correlation functions return "HBT radii" which fac- 
torize according to Equation[T6l i.e. F/^, which quantifies the dynamically-generated 
substructure, is identical in p + p and A+A collisions. However, the STAR data 
show significant non-femtoscopic structures which must be properly accounted 
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for (I32I) before drawing firm conclusions. If the factorization is unchanged after a 
more sophisticated treatment, the physics implications might be dramatic. We do 
not discuss this here, but simply observe that the NNUS scenario is a likely baseline 
expectation for p + p a\ 
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Fig. 18 Full markers: predicted pr-dependence of S,„,. and A fro m Gaussian fits to calculated 
correlation functions in a Pythia+hadronic rescattering model il60|). Here, the hadronization time 
was set to T = 0.1 fm/c, the value required to reproduce E735 jl6ll) data at the Tevatron. Also 
shown in open markers are the femtoscopic parameters coiTesponding to a longer hadronization 
time. In this case, there is little very hadronic rescattering and thus a much weaker -dependence 
of R„„,. 



The increase of HBT radii with multiplicit y has also been observed previously 
in p — p collisions by the E735 Collaboration ( Il6lh . While in A +A collisions, this 
is naturally related to incr easing length scales in the entrance channel geometry, 
Paic and Skowroriski (Il62l) postulate that jet dynamics, rather than bulk properties, 
drive this dependence in the p + p system. Within a simple model of hadronization, 
they can reproduce the E735 multiplicity dependence, and make predictions of sim- 
ilar multiplicity dependence for p + p collisions at the top LHC energy. However, 
since the expectation of increasing length scales with multiplicity seems to be rather 
generic to all scenarios, it will be interesting to see these predictions expanded to 
more differential measures- say, the multiplicity and pr dependence, probing both 
aspects of Equation [161 This should allow a more discriminating comparison be- 
tween models, allowing some to be ruled out. 

Such differential predictions have very recently been performed by Humanic (II60I) 
in the context of a Pythia-nhadronic rescattering (through HRM) scenario. It is found 
that, contrary to some expectations, hadronic rescattering is crucial to understand the 
M(8) pt dependence of the E735 data. Reproducing the data requires the assumption 
of a surprisingly short hadronization timescale (~ 0. 1 fm); longer timescales do not 
allow sufficient hadronic rescattering needed to describe the pj- -dependence. The 
prediction of the model for the highest multiplicity p + p collisions at ^ =1.4 TeV 
are shown in Figure [TS] 
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6 Summary 

From an extended and deformed initial state, to medium-induced energy loss of 
partons, to collective flow, to deconfinement itself, space-time considerations un- 
derlie any understanding of the bulk system created in a heavy ion collision. Trans- 
port model calculations that simulate the evolution of this highly dynamic system, 
produce space-time configurations which reflect basic physical quantities like the 
Equation of State, viscosity, and degree of equilibration. 

The most direct way to get at this space-time information is through femtoscopic 
measurements. With the availability of first nuclear beams a quarter century ago, 
a femtoscopic measurement was considered a success if it could demonstrate that 
the overall scale of the system was on the order of a few fermi. Since then, much 
improved experimental techniques and more precise data, coupled with advances in 
theory and phenomenology, have made femtoscopy a precision tool. A comprehen- 
sive program has resulted in multi-dimensional (e.g. "out-side-long") pion correla- 
tion functions as a function of y, pr and 0; similar measurements for other particle 
species combinations have been steadily accumulating. 

These systematic measurements allow the extraction of the size, shape, orienta- 
tion, timescales, and even (for non-identical particle correlations) the positions of 
the homogeneity regions. The momentum- and species-dependences of these quan- 
tities reveal the nature and magnitude of collective behavior through the signature 
geometric substructure it generates. 

Generically, the data indicate an emitting source whose overall space and time 
scales are determined by the final-state chemistry. A long-range component of the 
emission function has been identified and is consistent with expections from con- 
tributions from resonance decay. The strong space-momentum correlations probe 
longitudinal, radial, and elliptical collective flow. They reveal a system that is ap- 
proximately boost-invariant near midrapidity, dominated by a position-dependent 
collective boost field. For non-central collisions, the system's shape evolves away 
from its initial out-of -plane extended geometry, due to anisotropic flow fields cou- 
pled with finite evolution time. This flow field is the main feature reflecting the 
response of the bulk system to the extreme conditions of the collision. The flow- 
dominated scenario is quantitatively consistent with momentum-only measurements 
like pt spectra, rapidity distributions, and emission anisotropy (e.g. V2). However, 
any inference of the structure of the flow field, based on the momentum-only data 
alone, would be a faith-based exercise. 

Moving beyond the generic scenario requires comparison with transport calcu- 
lations. After initial consternation followed by detailed study, it now appears that 
the femtoscopic data set for heavy ion collisions from RHIC is on the verge of be- 
ing satisfactorily explained. The source of much of the mystery in understanding 
source sizes at RHIC derived from a conspiracy of model features, which had been 
neglected or poorly implemented in previous dynamic treatments. Although none 
of these were individually sufficient to explain the HBT puzzle, the conspiracy of 
several corrections to push predictions in the same direction appears to bring predic- 
tions and data into reasonable agreement. Morale and optimism are high, but signifi- 
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cant work remains to be done as the success has only been demonstrated for nn cor- 
relations. Simultaneous analyses of azimuthally sensitve HBT and of non-identical 
particle correlations are needed and underway. Furthermore, while toy parameteri- 
zations can simultaneously reproduce essentially all soft-sector data, it has yet to be 
demonstrated that femtoscopic observables can be explained with the same model 
as elliptic flow. Only then, can victory be declared. We are now, however, close. 

Transport models will be tested even more strenuously, as heavy ion data at the 
LHC becomes available. Several groups have generated femtoscopic predictions 
from Boltzmann and hydrodynamic calculations; thus far, they all pertain to two- 
pion interferometry only. Compared to those at RHIC, geometric scales are expected 
to increase on the order of 35%, though there are significant differences between the 
model predictions. All predict a break with the "universal" ^j- -dependence observed 
at RHIC and lower energies. In particular, they predict steeper fall-off of the trans- 
verse radii, due to much larger predicted transverse flow. This increased flow is 
further reflected in predictions of less-steep single particle spectra. 

The confidence with which we now state our understanding of heavy ion colli- 
sions may be threatened by similar data for small sources. Most of the same features 
and trends observed in heavy ion collisions, and often attributed to collective phe- 
nomena, persist for the smallest p+ p source sizes. An analogous situation may 
also exist for other observables: there are suggestions that, modulo finite number 
effects, chemical yields and species-dependent pj spectra from hadronic collisions 
are quite similar to those from heavy ion reactions. This challenges our understand- 
ing at both the quantitative and qualitative level, and might signal the need for a new 
theoretical paradigm for dynamics and hadronization of small systems. Certainly, it 
raises issues for further research at both RHIC and the LHC, where apples:apples 
comparisons are possible. 
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